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Topological insulators are a broad class of unconven- 
tional materials that are insulating in the interior but con- 
duct along the edges. This edge transport is topologically 
protected and dissipationless. Until recently, all existing 
topological insulators, known as quantum Hall states, vi- 
olated time-reversal symmetry. However, the discovery of 
the quantum spin Hall effect demonstrated the existence 
of novel topological states not rooted in time-reversal vi- 
olations. Here, we lay out an experiment to realize time- 
reversal topological insulators in ultra-cold atomic gases 
subjected to synthetic gauge fields in the near-field of an 
atom-chip. In particular, we introduce a feasible scheme 
to engineer sharp boundaries where the "edge states" are 
localized. Besides, this multi-band system has a large pa- 
rameter space exhibiting a variety of quantum phase tran- 
sitions between topological and normal insulating phases. 
Due to their unprecedented controllability, cold-atom sys- 
tems are ideally suited to realize topological states of mat- 
ter and drive the development of topological quantum 
computing. 



I. INTRODUCTION 

The search of exotic states of matter has been one of the 
key driving elements of condensed-matter physics. Among 
these unconventional phases, the quantum Hall effect has been 
the paradigm of topological phases not captured by Landau 
symmetry breaking. Recently, the subject of topological in- 
sulators has emerged as a new frontier, discovering novelties 
in the single-particle band theory and providing a new impe- 
tus to the constantly evolving quantum many -body physics of 
strongly correlated systems ^ The topological classification 
of insulators and superconductors points towards an inherent 
universality underlying these exotic states of mattei^. 

A new kind of topological insulator, based on the quantum 
spin Hall (QSH) effect^, was first theoretically predicted^ ^, 
then experimentally confirmed in HgTe/CdTe quantum wells^ 
and three-dimensional Bi2Te3^. Unlike quantum Hall effects, 
both integei^ (IQHE) and fractional^ \ the QSH effect occurs 
at zero magnetic field and thus respects time-reversal (TR) 
symmetry. The simplest model of the QSH effect consists of 
non-interacting charged spin- 1/2 fermions where the two spin 



components are described as IQHE states at equal but opposite 
"magnetic fields." This composite system produces degener- 
ate (Kramers) pairs of counter-propagating chiral edge states 
with opposite spins. Residing in the bulk gap of the band 
structure, edge modes are gapless excitations that provide the 
backbone for the topological classification of the insulators. 
In contrast to the usual IQHE where edge states carry electri- 
cal current with quantized conductance^^, the edge states in 
QSH systems carry robust spin currents. In the ideal model, 
the latter are associated to a quantized spin-conductivit}^^. 



In electronic materials QSH effects result from spin-orbit 
(S-0) coupling, and spin-mixing interactions destroy the 
quantization of the spin Hall conductivity. A simple way to 
determine the fate of QSH insulators upon spin mixing is to 
count the number of Kramers edge pairs inside the bulk gap. 
Due to TR symmetry, perturbations can only mix different 
Kramers pairs and gapless edge states are protected when their 
number is odd. This signals a QSH insulator whose properties 
are topologically different from a normal band-insulator (NI) 
that presents an even number of edge states-^^^. Therefore, 
a Z2 even/odd topological index captures the distinction be- 
tween QSH and NI quantum phases. The discovery of topo- 
logical QSH systems has opened the possibility to design new 
spintronic devices exploiting the spin-dependent direction of 
edge or surface states in order to manipulate and transport spin 
current without dissipation. 



Here, we propose an optical-lattice realizatiorPEIl of such 
quantum phases and show that cold atoms provide a simple 
system where quantum spin-Hall physics can be investigated. 
The origin of these QSH phases is rooted in engineered gauge 
rather than S-0 coupling. We describe an original 
scheme to synthesize gauge fields, exempt from spontaneous 
emission implicit in all previous proposals^^ The great ver- 
satility offered by cold atom experiments allows us to reach 
the purest realization of the QSH effect and to explore strik- 
ing aspects of this topological state of matter. Besides, the 
non-interacting limit - the ideal realization of QSH phases - 
and the stability of the topological phases can be explored by 
means of Feschbach resonances and engineered disordei^. 
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II. THE SYSTEM 



Realizing topological insulators with cold atoms is par- 
ticularly attractive i^^l^, and we hereby provide a concrete 
setup using fermionic ^Li subjected to synthetic gauge fields. 
There are numerous proposals for creating AbeliarPJES and 
non-Abelian^^ gauge fields, leading to interesting phenom- 
. The vast majority of these proposals depend on 
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laser-induced Raman coupling between internal atomic states; 
such a method was recently implemented with bosonic ^^Rb 
atoms creating an Abelian gauge field^^ For the alkali 
atoms, the requisite coupling limits the possible detuning of 
the Raman laser beams from the ground S to first excited P 
transition, about to the excited state fine- structure splitting 
[7.1 THz for ^^Rb, 1.7 THz for ^^k, and only 10 GHz for 
^Li]. The very small detunings possible for fermionic atoms - 
^^K and ^Li - imply large spontaneous emission rates and thus 
the desired gauge fields must be synthesized through an alter- 
nate scheme. Here we describe such a setup that combines 
state-independent optical potentials, with micron scale state- 
dependent magnetic potentials in an atom chip. Our discus- 
sion focuses on harmonically trapped quantum degenerate ^Li 
systems about 5 /im above the atom chip's surface. Figure [T] 
depicts how the required potentials for this proposal can be re- 
alized using a combination of static and radio-frequency mag- 
netic fields. By completely eliminating spontaneous emission, 
this approach makes practical the realization of gauge fields 
for all alkali atoms. 

Our proposal for realizing a fermionic model with a SU(2) 
gauge structure requires four atomic states \gi) = \F = 
1/2, m,. = 1/2), 1^2) = |3/2,-l/2), |ei) = |3/2, 1/2), 
and |e2) = |l/2, -1/2) [see Fig. [T] (a)], which are confined in 
a square optical lattice and described by the Hamiltonian 

^ = ~ t ^]n+l,n^ ^ C^rn,n+1^ ^C77^,n+h.C. 

m,n 

+ ^stag <^]n,n^m,n- (1) 

m,n 

Cm, n is a 2-component field operator defined on a lattice site 
{x = ma^ y = na), a is the lattice spacing, m, n are integers, 
and t is the nearest-neighbor hopping. In this model the ma- 
trices ^x,y modify the hopping along the x, ^/-directions and 
result from the synthetic gauge field [cf. below and Meth- 
ods]. Here, all the states experience a primary lattice potential 
Vi{x) = Vx sin^(A:x) along x which gives rise to a hopping 
amplitude t ^ 0.4 kHz. A secondary much weaker lattice 
V2{x) = 2Astag sin^(/cx/2) slightly staggers the primary lat- 
tice with Astag ^ t. These lattices, with an approximate period 
of 2 /im, are produced by two pairs of A = 1064 nm lasers, 
slightly detuned from each other and incident on the reflec- 
tive surface of the atom chip [see. Fig [T](b)]. Additionally, 
these beams create a lattice along z with a 0.55 /im period, 
confining the fermions to a 2D plane. 

We propose to design the SU(2) hopping operators as fol- 
lows 



where we have set the units a = h = 1, and <Jx,z are Pauli ma- 
trices. In order to engineer these state-dependent tunnelings, 
the states \g) and |e) must experience oppositely-signed lat- 
tices realized along the y direction of the chip, as in Fig.[T](c). 
This scheme exploits the Zeeman shift of an array of wires 
with alternating currents ±1. The Oy hopping operator in Eq. 
^ involves x-dependent phases that can be realized with an 
additional grid of wires spaced by a = 2 /im along x [cf. 
Methods] . Remarkably, these moving Zeeman lattices repro- 
duce the Raman-assisted hopping from previous schemes^^"^ 
and lead to the operator Oy. A potential gradient along y de- 
tunes this Raman coupling into resonance and is produced by 
simply shifting the center of the harmonic potential - equiva- 
lent to adding a linear gradient. The specific form of the hop- 
ping operator Oy is motivated by the fact that it corresponds 
to opposite "magnetic fluxes" ±a for each spin component, 
where a = p/q and q are mutually prime numbers. Finally, 
the additional contribution to the hopping that mixes the 
|e) and \g) states can be realized in a like manner. The cor- 
responding control parameter 7, and the additional staggered 
potential which induces alternate on-site energies e = ±Astag 
along the x-axis, are shown below to enrich the phase transi- 
tions in a novel manner. Let us emphasize that the Hamilto- 
nian ([l])-(|2]) satisfies Time-Reversal invariance, since it com- 
mutes with the TR-operator defined as T = id-yK, where 
K is the complex-conjugate operator. This fundamental sym- 
metry opens the possibility to realize the first instance of Z2- 
topological insulators in a cold-atom laboratory. Finally, we 
stress that the two components of the field operators corre- 
spond in general to a pseudo-spin 1/2, but in our proposal refer 
specifically to spin components of ^Li in its electronic ground 
state. 

In this paper, we present results for the case a = 1/6, 
namely when the synthetic gauge field ^ corresponds to 
Oy{m) = 27rmaz/6. This choice exhibits an extremely rich 
phase diagram with almost all possible topological phase tran- 
sitions [cf. Figs. |4|. The topological phases discussed be- 
low are robust against small variations Sa ^ 0.01: indeed the 
properties rely on the existence of bulk gaps which are contin- 
uously deformed when a is varied, as in the Hofstadter butter- 
fl>^^. Moreover, other configurations of the gauge field could 
be experimentally designed and would lead to similar effects. 



In Sect. Ill we investigate the phase transitions between topo- 
logical insulators for a cylinder geometry in the absence of the 



27r7ax, 



and 



Oy = 27rxaaz, (2) 



confining potential Konf- We then show in Sect. [IV| how these 
properties can be detected when this confinement is applied. 



III. TOPOLOGICAL INSULATORS AND PHASE 
TRANSITIONS 

When 7 = 0, Eqs. ([l])-(|2]) describe two uncoupled quan- 
tum Hall systems. For generic a = p/q, the fermion band- 
structure splits into q subbands, the lattice analogues of Lan- 
dau levels. These subbands are separated by gaps, or by 
"pseudogaps" exhibiting nodal Dirac fermions. Our setup 
thus provides a multi-band system, where a variety of band- 
insulators, separated by metal or semi-metal regions, can be 
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Figure 1: Experimental realization using an atom-chip and resulting model (a)-(d): (a) Level diagram identifying the four required states 
from within the ground 2Si/2 electronic manifold of fermionic ^Li. Each of the state-pairs {l^'i), \g2)} and {|ei), |e2)} are clock state pairs 
(with the same magnetic moment) and the indices 1 and 2 label the SU(2) degree of freedom. Moreover the \gi,2) and |ei,2) pairs have 
opposite magnetic moments, so they experience opposite-signed magnetic potentials, (b) Explicit experimental layout showing the origin of 
optical (state independent) and magnetic (state dependent) potentials and coupling fields. First, a state-independent, staggered, lattice along x is 
formed by the separate interference of two pairs of A = 1064 nm laser beams slightly detuned from each other to eliminate cross interference. 
The respective intersection angles are chosen so the lattice period differ by a factor of two. Both beams reflect from the chip-surface and 
form vertically aligned lattices, trapping the degenerate Fermi gas about 5 /im above the surface. The inset shows the chip geometry, from 
top to bottom: a reflective chip-surface, gold wires aligned along y with a 1 /xm spacing along x (producing the \g)-\e) coupling), and finally 
gold wires aligned along x with a 2 /xm spacing (producing the state dependent lattice), (c) Computed atomic potentials and coupling matrix 
elements, (d) The resulting lattice with couplings schematically illustrated. 



reached by varying the atomic filling factor. As discussed 
below, some of these insulators are topologically non-trivial 
and feature gapless edge states. The latter are localized at the 
boundaries of the sample, and correspond to gapless excita- 
tions. When the Fermi energy Ep lies inside a bulk gap, the 
presence of these states is responsible for the spin transport 
along the edges. 

Figure [2] (a) illustrates the energy spectrum obtained for a 
cylindrical geometry [cf. Methods]. In this figure, the bulk 
bands [thick blue bands] are clearly differentiated from the 
gapless edge states within the bulk gaps [thin purple lines]. 
In the lowest bulk gap, the edge state channel A, and its TR- 
conjugate B, correspond to localized excitations which travel 
in opposite directions [cf. Fig|2](b)]. Each boundary is popu- 



lated by a single pair of counter-propagating states, each cor- 
responding to an opposite spin [i.e. helical edge states]. Ac- 
cordingly, this lowest bulk gap describes a topological QSH 
phase. Conversely, the next gap located at E ^ —1 is tra- 
versed by an even number of Kramers pairs, thus yielding a 
normal band insulator. Remarkably enough, different topo- 
logical insulators can be experimentally engineered in this 
multi-band scenario by simply varying the atomic filling fac- 
tor. The Fermi surface at half-filling, namely Ep = 0, corre- 
sponds to a discrete set of isolated points, and thus describes 
a semi-metal phase. 

An alternative approach to the above even-odd criteria re- 
lies on the computation of a Z2 topological invariant that 
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Figure 2: Bulk energy bands and edge-state channels: (a) Energy 
spectrum of the uncoupled system (7 = and a — 1/6) computed 
in a cylindrical geometry: the bulk energy bands [thick blue bands] 
are traversed by edge states [thin purple lines]. The energy is in units 
of the tunneling amplitude t. (b) Schematic representation of the 
two edge states, A and B, that lie inside the first bulk energy gap 
depicted in (a). The spins traveling around the edges are respectively 
represented by red and green arrows. 



characterizes the bulk gaps [see IVEethods]. The QSH phase 
is distinguished from the normal insulator by the Z2 index 
■> which respectively takes the value Iz^ (QSH) = 1 and 
7^2 (NI) = 0. We computed this index inside the four gaps 
depicted in Fig|2](a), obtaining the following sequence 

Jz, ={1,0,0,1}, 

in perfect agreement with the above edge- state analysis. 
Therefore, in the uncoupled case 7 = where spin is a 
good quantum number, the spin conductivity is quantized as 
cTg = 2 ^ is related to the spin current — Jt ~ h- 
Additionally, this quantized spin Hall conductivity is equal to 
the difference of the Chern numbers associated to the up- and 
down-spin [cf. IVLethods]. 

In our multi-band system, the lattice-potential distortions 
can drive direct transitions between normal and topological 
insulating states. Figure [3] shows the bulk gaps and edge states 
for increasing values of the staggered potential's strength 
Astag. In the simple case 7 = 0, where two quantum Hall sys- 
tems are superimposed, we already obtain various realizations 
of QSH phases. By tuning Astag, one can observe a NI to semi- 
metal transition, and more importantly, a semi-metal to QSH 
phase transition in the bulk gaps located around E = ±.1. 
Remarkably enough, we have shown that the staggered poten- 
tial, which can be controlled in an optical-lattice experiment, 
induces a quantum phase transition from a NI to a QSH phase 
at the critical value Astag = ^ [cf. Fig.[3](b)]. At half-filling, 
the existence of Dirac points appear to resist the opening of 
the gap for small staggered potential and eventually lead to a 
NI phase for Astag > 1.25t [cf. Fig. [3](d)]. As we describe 
below, the latter situation radically changes as the coupling 7 
is switched on. 
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Figure 3: Energy bands E{ky) for 7 = and a = 1/6 
, in the presence of an external staggered potential (a)-(d): (a) 

Astag = 0.5t, (b) Astag = t, (c) Astag = 1.25t, (d) Astag = 1.5t. The 
arrows point at the open bulk gaps and indicate their corresponding 
topological phases. The purple rectangles highlight the NI to QSH 
phase transition that occurs around Ep — 1. The energy is expressed 
in units of the tunneling amplitude t. 

The parameter 7 7^ perturbs the ideal quantum spin Hall 
system and fundamentally changes the aforementioned phase 
transitions. The Z2 invariant analysis provides a particularly 
useful and efficient tool to obtain the full phase diagram of 
the system in the (7, Astag) parameters plane [cf. Methods]. 
The phase diagram represented in Fig.|4ja) has been obtained 
numerically by evaluating the Z2 index inside the bulk dXE ^ 
±1 in the vicinity of 7 = 0. Here we observe three distinct 
regions, namely the metallic (blue), QSH (red), and NI (cyan) 
phases. These three phases coexist at a tricritical point situated 
at 7 = and Astag = ^- It is interesting to highlight that the 
QSH phase occurs for a wide range of the parameter 7, which 
indicates the robustness of this topological phase under small 
local perturbations of the Hamiltonian ([T]). 

The wide range of possibilities offered by optical-lattice ex- 
periments enables us to consider the strong coupling regime 
corresponding to 7 = 0.25. In this limit, the previously in- 
dependent quantum Hall subsystems (7 = 0) become maxi- 
mally coupled, a fact that drastically modifies the topological 
phase transitions presented above. In Fig. [5] we illustrate the 
bulk bands and edge states for increasing values of the stag- 
gered potential. Each gap is associated to a QSH or a NI phase 
according to the even-odd number of TR pairs. In this strong- 
coupled regime, we observe a completely different phase tran- 
sition associated to the bulk gap at half filling. Contrary to the 
semi-metal to NI phase transition described in the uncoupled 
case [cf. Figs.[3ja)-(d)], we identify here a QSH to NI tran- 
sition which occurs after a gap-closing at Astag = 1.25t. This 
process is emphasized on the figures with purple circles [cf. 
Figs.[5](b)-(d)]. The QSH phase can only occur at half-filling 
for large couplings 7. Furthermore, one also notices that the 
neighboring bulk gaps at £^ ^ ±1 present exactly the oppo- 
site phase transition, from a NI to a QSH phase. The latter 
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show that the parameters 7 and Astag play competing roles: 
7 tends to shrink the gaps width, while Astag tends to open 
them. In this multi-band configuration, an opening of one of 
the gaps may result in a gap-closing of the neighboring bands. 
Therefore, by manipulating the coupling between the states, 
the staggered potential and the Fermi energy, the system dis- 
cussed here offers the possibility to explore different topologi- 
cal phase transitions within the several bulk gaps. These novel 
features endow the topological phase diagram with an intrin- 
sic richness and complexity, not present in other condensed- 
matter realizations of the QSH effect. 
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Figure 4: Phase diagrams in the (7, Astag) plane (a)-(c): (a) In the 

vicinity of the uncoupled case 7 = when Ef = ±1. The tricritical 
point tc is situated at (7 = 0, Astag = t) and corresponds to Fig. [s] 
(b). (b) In the vicinity of the maximally coupled case 7 = 0.25, 
when Ef — 0: two tricritical points ti and t2 are observed, (c) In the 
vicinity of the maximally coupled case 7 = 0.25, when Ef — ±1- 
The red, blue and cyan regions respectively show the QSH, metallic 
and normal insulating phases. The staggered potential's strength Astag 
is expressed in units of the tunneling amplitude t. 



occurs at A. 



t, as emphasized on the successive figures 
with purple rectangles [cf. Figs.[5](a)-(c)]. Both Figs. |3] and |5] 



Figure 5: Energy bands E{ky) for 7 = 0.25 and a = 1/6, 
in the presence of an external staggered potential (a)-(d): (a) 

Astag = 0.5t, (b) Astag = t, (c) Astag = 1.25t, (d) Astag = l-5t. The 
arrows point at the open bulk gaps and indicate their corresponding 
topological phases. The purple rectangles [resp. circles] highlight 
the NI to QSH [resp. QSH to NI] phase transitions. The energy is 
expressed in units of the tunneling amplitude t. 
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In our model, topological phase transitions of different na- 
ture can exist in each of the bulk gaps. As shown in Fig. |5] 
neighboring bulk gaps can undergo reverse transitions succes- 
sively. To fully capture the richness of this phenomenon, we 
numerically compute the Z2 index for a wide range of the pa- 
rameters around 7 = 0.25. Figure |4jb)-(c) show the bulk gaps 
at £^ = and E ^ ±.1, respectively. In Fig.|4jb), the central 
gap features tri-critical points, and the QSH-NI phase transi- 
tions occur along a well-defined curve. On the other hand, 
in the neighboring gap [Fig.Qc)], the QSH to NI phase tran- 
sition is separated by an intermediate metallic region. There- 
fore, we conclude that a single realization of this model allows 
us to study a wide variety of different topological phase tran- 
sitions. 
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IV. DETECTION IN THE PRESENCE OF THE 
CONFINING POTENTIAL 

A. The harmonic trap and gauge fields with anisotropic 
hopping 



We now describe a new feasible scheme to engineer a sharp 
interface where edge states can be locaHzed. This is essential 
for detecting topological states in optical lattices, where the 
indispensable harmonic trap used to confine atoms strongly 
influences the edge states^^ when Konf (edge) ~ A, where 
A is the gap's width. In the present context A ~ 0.5 — It, 
requiring temperatures T < 10 nK. 



The key aspect of our proposal exploits the fact that the hop- 
ping along the ^/-direction, ty, is laser-controlled and hence 
can be tuned. Furthermore, since the harmonic trap has a min- 
imal effect at the center of the trap, we divide the chip into 
three regions: the central region is characterized by a large 
hopping ty, while the two surrounding regions feature small 
hoppings ty <^ ty. This highly anisotropic hopping creates a 
sharp interface where the edge states of the central region - a 
topological insulator - localize. In other words, by controlling 
the strength of the Raman coupling, we squeeze the energy 
bands describing the outer parts so that it does not interfere 
with the bulk gaps of the central topological phase where the 
edge states reside. We show that the helical edge states are 
robust and exactly localized at the interface between the inner 
and outer regions, whose location can be easily controlled in 
this scheme [cf. Fig.[6](a)]. Furthermore, since the topological 
phases are confined in the center of the trap, one verifies that 
the phase diagrams discussed in the previous sections are valid 
for a much wider range of the harmonic potential's strength. 
These observations are illustrated in Fig. [6j where the energy 
spectrum and edge- state wave-functions are depicted for dif- 
ferent configurations. In Fig.[6](b), we observe a new metal- 
lic continuum band corresponding to the outer region, which 
corrupts all but the lowest and highest band gaps. When the 
Fermi energy lies in the preserved gap, we see how the edge 
states are robustly localized within the interface [cf. Fig. |6] 
(a)]. Moreover, Fig. |6](c) shows that the staggered potential 
opens a bulk gap at half-filling, allowing to explore the phase 
diagrams of Fig.|4](b)-(c). 



B. Detecting the edge states and the quantum spin Hall effect 

Probing the topological edge states in a cold-atom experi- 
ment remains a fundamental challenge^^^^"^^. Although den- 
sity measurements might already detect the accumulation of 
atoms at the designed edges [i.e., the aforementioned in- 
terfaces], alternate techniques should access more profound 
properties associated to these specific states. For instance, 
stimulated Raman spectroscopy techniques can measure the 
dispersion relation of these excitations^^. A further promising 
possibility would be to exploit light Bragg scattering in or- 
der to observe signatures in the dynamical structure factoi^^ : 
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Figure 6: Edge-states for the anisotropic-hopping scheme in the 
presence of the trap (a)-(c): (a) The chip is divided into three 
regions: two sharp interfaces separate the inner region with hop- 
ping amplitude t'y from the outer regions [with hopping amplitude 
ty ^ t'y\- The wave-functions amplitudes are localized 

at the interfaces and correspond to the eigenstates dX E — [cf. 
Fig.[6](b)]. The red [resp. green] amplitudes correspond to up- [resp. 
down-] spins, and they are filled [resp. empty] according to their 
orientation along the "edges". This color code highlights the helical 
nature of the "edge-states" at the interfaces, (b) Energy spectrum 
E — E(ky) for the case 7 = Astag = in the presence of a har- 
monic potential, K:onf(edge) = 0.5t. The hopping parameters of the 
inner and outer regions are respectively t'y — 2t and ty = O.lt. The 
eigenvalues are colored in black, blue or purple when their corre- 
sponding eigenstates are localized at the interface, the inner region 
or the outer regions respectively, (c) Energy spectrum E — E{ky) 
for the case 7 = 0.25, Astag = 1.125t in the presence of a har- 
monic potential, K:onf(edge) = 0.5t. The hopping parameters of the 
inner and outer regions are respectively t'y — lAt and ty — O.lt. 
Note how the edge- states channels are robust when the anisotropic- 
hopping scheme is applied to the trapped system. 
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whenever the Fermi energy crosses a gapless edge-state, a 
delta peak is distinctively observed in the structure factoi"^^. 
The present proposal is particular suited for testing the latter 
detection method, since the light scattering can be focused on 
the interfaces where the edge- states have been shown to be 
robust in the presence of a confining harmonic potential. 

Finally, supplementary signatures of the quantum spin Hall 
effect could be detected through spin-resolved density mea- 
surements^^. A synthetic electric field [e.g. by accelerating 
the lattice or using gravity], should separate the different spin- 
components, thus detecting spin-currents in the transverse di- 
rection. In the 7 = case, the spin-Chern number is given 
by SChN=A^| - = [see Methods], which can be 
evaluated through the Streda formula applied to the spin-up 
densit}!^. 



V. CONCLUSIONS AND OUTLOOK 

In this article, we described a concrete and promising pro- 
posal of synthetic gauge fields in square optical lattices that 
overcomes the severe drawbacks affecting earlier schemes. 
Moreover, we showed how such engineered gauge fields are 
ideally suited to experimentally realize the most transpar- 
ent QSH phase, simplifying the usual S-0 induced effects in 
condensed-matter systems. Accordingly, ultra-cold atoms of- 
fer a rich playground to study time-reversal topological states 
of matter. 

In systems where SU(2) gauge fields encode a non-trivial 
topology in the band structure, we showed how a staggered 
potential can drive a wide range of quantum phase transitions 
between metalHc, QSH and NI phases. Both NI-QSH and 
QSH-NI transitions occur successively in different bulk gaps 
for different values of the staggered potential. Such multi- 
band gap-dependent quantum phase transitions constitute a 
unique and rich feature of our system. 

The cold-atom realization of topological band-insulators 
and helical metals proposed in this paper will pave the way 
for engineering correlated topological superfluids and insu- 
lators. Superfluids arising from pairing instabilities in he- 
lical metals are of particular importance for quantum com- 
puting due to the existence of Majorana fermions in vor- 
tex cores^^^^. Highly controllable cold atoms are not only 
better suited than condensed-matter systems for manipulat- 
ing vortices and quantum information stored in their Majo- 
rana modes, but also have slow dynamics, making them very 
promising for quantum computing applications. 

Finally, let us comment on the possibility of synthesizing 
SU(N) gauge fields with N > 2. Considering fermionic atoms 
with higher hyperfine levels F > ^, it is possible to create 
artificial gauge fields where N = 2F + 1 is even. In this re- 
spect, one can envisage situations where the helical edge states 
present a richer spin structure, and thus offer the opportunity 
to explore new avenues and exotic topological phases. 



VI. METHODS 

A. Synthetic gauge fields and the quantum spin Hall effect 

In condensed-matter systems, the quantum spin Hall ef- 
fect finds its origin in the inherent spin-orbit coupling of 
the materials^ ^. In the seminal lattice model of Ref."^, this 
term involves next-nearest-neighbor (NNN) hoppings on a 
honeycomb-lattice 

%o = ^Aso ^ lyijclcTzCj, (3) 

{{iJ)) 

where = ±1 according to the orientation of the path sepa- 
rating the NNN i and j. In analogy to the Haldane model for 
spinless electrons^^, this term can be expressed as 

Hso = Aso E ^e^^^'-'^c,, (4) 

where A. = Aaz is a gauge field defined on the two trian- 
gular sub-lattices, and such that A. - d\ = ±f accord- 
ing to the orientation of the path. Note that on a lattice, the 
gauge field A. modifies the hopping of a particle through the 
standard Peierls substitution^^. In the original model of Hal- 
dane [i.e. v4 = A and A d\ = ±f ], this term opens 
an integer quantum Hall gap, while in the spin- 1/2 model of 
Ref."^, this term is proportional to the Pauli matrix and thus 
leads to a quantum spin Hall gap. Similarly, in the GaAs semi- 
conductor system of Ref.^, the spin-orbit coupling is described 
by a gauge field A = A(Jz, where V x A is analogous to a 
uniform magnetic field. In both models, the presence of the 
Pauli matrix in the gauge field A highlights the fact that 
quantum spin Hall systems can be interpreted as the coupling 
of two quantum Hall subsystems. 

In the present setup, one proposes to engineer the coupling 
of two Hofstadter subsystems. The Hofstadter system de- 
scribes spinless electrons on a square lattice subjected to a 
uniform magnetic field and leads to a multi-gap energy spec- 
trum characterized by non-trivial topological invariants^. Ac- 
cordingly, our square lattice Hamiltonian ([T]) features the hop- 
ping matrices ^x,y, which are controlled in our optical-lattice 
setup, and which can be associated to the presence of a syn- 
thetic gauge field A as follows 



<9x= / A-dl, Oy= A'dl, (5) 

where we have set = c = e = 1. Here, the synthetic 
gauge field A is chosen in order to induce the ^x,y operators 
introduced in Eq. (|2]). The specific form of the hopping oper- 
ator Oy (X adz corresponds to opposite "magnetic fluxes" ±a 
for each spin component and thus leads to quantum spin Hall 
gaps. 
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B. Experimental proposal 

Here we describe our technique to create an atomic en- 
semble subject to a SU(2) gauge field (focusing on ^Li, but 
equally applicable to ^^K) that combines state-independent 
optical potentials, with micron scale state-dependent magnetic 
potentials near the surface of an atom chip. Our discussion 
focuses on harmonically trapped quantum degenerate ^Li sys- 
tems a distance h = b /im above an atom chips surface, in 
the presence of lattices with period a = 2 /im along x and y. 
We use hkL = irh/a and El = h^k\/2m = h x 8.3 kHz - 
the recoil momentum and energy momentum - as the natural 
dimensions of energy and momentum. 

All proposals to create light-induced gauge fields for ultra- 
cold neutral atoms somehow involve Raman transitions that 
couple ground states of unlike spin that differ in momentum 
by hSk by way of intermediary (and in practice off-resonant) 
excited states. For a typical system with two ground states, the 
resulting effective Hamiltonian has (potentially many) terms 
of the form 

, 1^ / exp[z(Jk-x + (/))]\ 

^ ^""^ ~ 2 V exp [-i{Sk • X + ^)] J 

= ^ [cos ((5k • X + ^) ax — sin {Sk • x + ^) ay] . 

Formally this is the same as the linear Zeeman shift gfl- B for 
a magnetic field in the x-y plane orbiting about z with spatial 
period 27r/|^k|. We take advantage of this formal similarity 
and constructs the requisite magnetic field near the surface of 
an atom chip. 

Our proposal for realizing an atomic model with a SU(2) 
gauge structure requires four atomic states with very specific 
properties outlined below. These physically correspond to hy- 
perfine states, labeled by |F, m/), in the ^Li 2S electronic 
ground state: \gi) = |l/2, 1/2),|^2) = |3/2, -l/2),|ei) = 
|3/2, 1/2), and les) = |l/2,-l/2). For simplicity, the la- 
bels 1^) and |e), have been preserved from Ref.'^ where they 
denoted physical ground and excited states. 

(A) The atomic states must all experience a primary lat- 
tice potential Vi{x) = VxSm'^{kx) along x; a lattice depth 
Vx = QEl ^ hx bO kHz produces a hopping matrix element 
t ^ O.Ob El ^ h X 400 Hz. A secondary, much weaker lat- 
tice V2{x) = 2Astag sin^(/cx/2) slightly staggers the primary 
lattice, where Astag ~ t <C Kc. Figure [l](b) shows two pairs 
of A = 1064 nm lasers, slightly detunea from each other, in- 
cident on the atom chip's reflective surface, producing lattices 
with period A/2 sin ^. For our base a = 2 fim lattice period, 
these angles [shown in Fig. [T| are Oi ^ 15° and 62 ~ 7.5°. 
These beams also produce a lattice along z with a 550 nm 
period, confining the fermions to a 2D plane. 

(B) The states \g) and |e) experience oppositely- signed lat- 
tice potentials along y as in Fig. [T] (c). This can be directly 
implemented with the Zeeman shift ^/i^ |^| of atoms pro- 
vided that 1^1,2) and |ei,2) have equal, but opposite magnetic 
moments ^ (i.e., are each on clock transitions). For the states 
- 1^1), 1^2), and |e2) - shown in Fig.[T](a) the magnetic 
moments are correctly signed and differ by less than 1% in 
magnitude at a bias field B = 0.2b G. 



With these states, a state-dependent lattice potential can be 
generated by an array of current carrying wires with alternat- 
ing +/ and — / currents, spaced by a distance a (repeating 
with full period 2a = 4 /im); because the evanescent mag- 
netic field decays exponentially ex exp (— 27r/i/a), the height 
h of the atoms above the surface must be comparable to a. A 
modest I = b /iA current^"^ in wires 3 /im below the chip- 
surface (8 /im from the atoms) produces slQ El Zeeman lat- 
tice, with a negligible 3 Hz hopping matrix element. 

(C) Our model also requires hopping along y with an x 
dependent phase: tyex.p{iqx). This can be realized with 
an additional grid of wires spaced by a = 2 /im along x, 
with currents Im- This provides moving Zeeman lattices with 
wave-vector q (leading to effective "Raman coupling") when 
In = Iosm{qma - ut). The a;/27r ^ 228 ± 0.23 MHz 
transitions indicated with arrows in Fig. [T] (c), |ei) \gi) 
and 1 62) 1^2), are independently controllable in phase, am- 
plitude, and wave-vector by commanding concurrent running 
waves at the indicated resonant frequencies. For this study 
we confine ourselves to the equal amplitude and wave-vector 
case. The minimum wavelength d = 27T/q of this moving 
lattice is Nyquist limited by d > 2a. 

In the frame rotating at the angular frequency uj, and after 
making the rotating wave approximation the coupling terms 
have the desired form tex.p{iqma). In the model Hamilto- 
nian ([T]), the phase for hopping along ^ is a = qa /27r = a/d 
in terms of physical parameters. In this proposal we focused 
on a = 1/6, implying d = 6a far exceeding the 2a Nyquist 
limit. 

Our scheme also requires, a contribution to the hopping 
along X that mixes the |e) and \g). Like in (C) this can be 
realized using a Zeeman lattice moving along x, but tuned to 
drive transitions between \gi) \g2) and |ei) {62). 

(D) To create effective edges within the trapped cloud of 
atoms, we propose to abruptly change ty from a larger value in 
the bulk to a smaller value ty as a function of x. Here, creating 
such a change in ty can be simply implemented by making the 
currents Iq in our wires x-dependent such that they sharply 
change from a large value in the center to a smaller value near 
the systems edge, as is schematically shown in Fig. 6 (a). 

(E) Lastly, a potential gradient along y detunes into reso- 
nance the Raman fields providing the y tunneling and is pro- 
duced by simply shifting the center of the harmonic potential 
- equivalent to adding a linear gradient. This gradient implies 
that the coupling fields uj^ will be slightly different in fre- 
quency to be simultaneously resonant. We select a gradient 
so the site-to-site (spaced by a) detuning is large compared 
to both the ^ 3 Hz bare and ^ 400 Hz Raman stimulated 
hopping along y, but small compared to the ^ 20 kHz band 
spacing along y. For example a 1 kHz/a = 500 Hz//im gra- 
dient results from shifting a typical 40 Hz harmonic trap by 
500 /im. This can be effected more efficiently in the labora- 
tory by using the linear gradient at the side of gaussian laser 
beams. 
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C. The cylindrical geometry description 

In order to compute the edge states, we consider a cylin- 
drical geometry with boundaries at x = and x = L. 
In this case, we apply periodic boundary conditions along 
the ^/-direction, and consider the partial Fourier transform 
Cm,n = Xl/c e~^^y^Cky{m). The Hamiltonian ([T]) becomes 
a collection of independent tight-binding models along the x- 
axis H = Y,ky HjBiky), where 

HjBiky) =-t^Cky{m^ l)te^^-Cfc^(m) +h.c 

m 

-^^^ky{myD{m)cky{m), (6) 

m 

and where the diagonal matrix yields 

D{m) = — 2t diag(cos(27ram + ky) ^cos{27ram — ky)) 

+ Astag(-l)"^I. (7) 

The energy spectra computed with this method presents the 
bulk bands as well as the edge states lying inside the bulk 
energy gaps [cf. Fig.[2|. 

D. The Z2 index and the spin-Chern numbers 

In analogy to the topological characterization of the inte- 
ger quantum Hall effect by means of Chern number J^^''^, we 
evaluate the Z2 invariants of the QSH effect through the com- 
putation of the ''spin-Chern-number'\SChN) ^3 45 latter 
topological invariant is associated to the bulk and is there- 
fore defined on the torus obtained by closing the cylindrical 
geometry [i.e. by identifying the two edges at x = and 
X = L]. Remarkably, the non-triviality of the SChN witnesses 
the presence of gapless edge-states, hence building an elegant 
bridge between the open and closed geometry descriptions 

The computation of the SChN requir es the consideration 
of twisted-periodic-boundary conditions^^EUSl ^ particle 
crosses the boundary located at x = = L, it acquires an 
extra phase e*^" , where a =t, | is the spin index. Since the 
spectral bulk gaps and their correspondong topological order 
are immune to the boundary conditions, the topological SChN 
is obtained as the average of the Berry's curvature over all 
the phases aquired through the boundary. Here we set 0^ = 
—Oi = qO for convenience and the twist-space spanned 
by {0,ky) describes a torus: G [0, ^1 ^nd ky G [0,27r]. 



Note that these boundary conditions break the translational 
symmetry along the x-direction, while preserving it along the 
^/-direction. Besides, the closing of the cylinder is consistent 
as long as L = integer x q. The single-particle Schrodinger 
equation can then be expressed in the form of a Harper-like 
difference equation 

E 

Ri^m+l ^m-l ^D{m)ll)m = - y V^m , ^ 7^ 1 , ^, 

Re'^^^^^Ji + R^ ^L-i + D{L)^L = -J^L. (8) 

where iprn is a two-component spinor, R = (cos(27r7) I + 
i sin(27r7)<Ja^) and D{m) is defined in Eq. [t] 

The topological invariant associated to the r* bulk energy 
gap, SChN(r), is computed on a discrete Brillouin zone ac- 
cording to the expression 

SChN(r) = ^^^(*(^,A:,)). (9) 

Here J-^(^) represents the non-Abelian Berry's curvatur^ 
associated to the multiplet ^ = (/)2, . . . , ^at), where 
(/)A denotes the single-particle wave-functions associated to 
the energy bands Ex lying below the r* bulk gap. The NI 
phase [i.e. = 0] is reached in the bulk gap when 

SChN(r) = mod4, while this gap corresponds to the QSH 
phase [i.e. Iz^ = 1] when SChN(r) = ±2mod4. 

In the uncoupled regime 7 = 0, our setup displays a dou- 
ble quantum Hall system in which particles with opposite 
spins are coupled to opposite magnetic fields, but not to each 
other. Such a system exhibits TR-symmetry, and therefore, 
bands come in time-reversed pairs characterized by opposite 
Chern numbers: ChN^ = — ChN^ . The net Chern number 
ChN = ChN^ + ChN| which accounts for the charge trans- 
port is zero, but the difference SChN = ChN^ — ChN^ is 
finite and defines a quantized spin Hall conductivity. The suc- 
cessive energy bulk gaps shown in Fig. [2]are associated to the 
Chern numbers [from top to bottom] ChN^ = { — 1,— 2,2,1} 
in agreement with the computed Z2 indices = {1, 0, 0, !}• 
However, bulk gaps associated to even Chern numbers lead to 
a vanishing spin Hall conductivity, as such a system has even 
pairs of edge states. Therefore, in systems with TR symmetry, 
Z2 invariants are the key to distinguish QSH and NI phases. 
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